Axially-symmetric Neutron stars: Implication of rapid rotation 
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Models of relativistic rotating neutron star composed of hyperon rich matter is constructed in the 
framework of an effective field theory in the mean-field approach. The gross properties of compact 
star is calculated at both static and the mass-shedding limit in the axially symmetric basis. The 
effect of appearance and abundance of hyperons on equation of state of dense matter and stellar 
properties is lineated with particular emphasis on the underlying nuclear interactions. We find 
that the models can explain fast rotations, which supports the existence of millisecond pulsars. An 
important offshoot of the present investigation is that, irrespective of the model parameters and 
interaction taken, the star seems to sustain faster rotations (an increase in rotational frequency up 
to « 50%) without any further deformation. 
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I. INTRODUCTION 

Neutron stars are the manifestation of the densest form 
of cold matter in the observable universe [l|, Q • Formed 
in the aftermath of type II supernovae, they represent 
the end point of the life cycle of the star. Since all the 
known forces of nature amalgamate to determine the evo- 
lution and structural properties of these stars, they can 
provide vital link between astrophysics, nuclear physics, 
particle physics as well as heavy-ion collision data. Theo- 
retically, the properties of neutron stars such as its mass, 
radius and related properties are the imprints of a par- 
ticular equation of state (EoS). At « (2 — 3)poj where po 
is the normal nuclear matter density, the presence of ex- 
otic forms of matter like hyperons, quark or mixed phase 
starts to appear and has substantial impact on the EoS, 
which in turn is reflected in the global properties of the 
star. Since we are still devoid of having any smoking gun 
signals from observations pertaining to the presence of 
these novel phase, apparently it becomes important to 
investigate and address these issues in a proper theoret- 
ical framework. It is known that the maximum mass of 
the star is regulated by the EoS well beyond 2po and the 
radius is regulated by the EoS in the immediate vicinity 
of pq. Arguably, in order to minimize the uncertainty in 
the mass and radius predictions, one needs to constrain 
the EoS both at po and beyond. Till date, more than 
1800 pulsars have been discovered with varying degrees 
of predictions, which invoke both challenge as well as in- 
terest in search for the realistic EoS for neutron stars. 
For example, on one hand, the evidence from the white 
dwarf-pulsar system JG751+1807 with a mass of M 
(2.1 ± O.2)M0 (Icr) and M = (1.6-2. 5)Mq {2a) [3] and 
a mass limit M> (2.10 ± O.28)Af0 from x-ray bursts of 
the neutron star EXO 0748-676 [3] rule out the possibil- 



ity of a soft EoS. On the other hand, the bounds from 
thermally emitting neutron stars eg. RXJ 1856.5-3754 
^ imply the star radius to be beyond 10-12 km empha- 
size the need of a stiff EoS. Similarly, the subthreshold 
production of kaons as measured by the KaoS collabora- 
tion \^ were investigated in transport models [t} which 
predicts the EoS to be soft. Also, recent results obtained 
from transport models reveals a soft symmetry energy at 
Pb > 2po [Sl has raised crucial implications on the na- 
ture of EoS. These apparent contradictions leaves ample 
opportunity to probe the various aspects of an EoS in 
detail and emphasize the need of constraints. 

With this motivation, we now employ models of rela- 
tivistic effective field theory Q to address these inevitable 
issues relating to the neutron star properties and dis- 
cuss implications of fast rotations on the stellar structure. 
The model embodies various interactions and cross cou- 
plings, that satisfies the nuclear matter saturation prop- 
erties and exploits the poorly known density dependence 
of the symmetry energy [lO, [IH. The model is gener- 
alized to include the lowest lying octet of baryons and 
is employed to study the static as well as rotational at- 
tributes of the star. 

The paper is organized as follows: In section II we 
give a brief description of the hadronic model which we 
use in our present work. In section III, we describe the 
numerical scheme to calculate the global properties of 
compact star. The main features of EoS and the global 
properties of the rotating neutron star is presented in 
section IV. We compare our results with other neutron 
star models in the light of observational bounds in section 
V. Finally, we conclude with the important findings of the 
present work in section VI. 



II. THE EQUATION OF STATE 
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In the framework of Quantum hadrodynamics (QHD) 
for describing the nuclear many-body problem, as a rela- 
tivistic systems of baryons and mesons, we employ an ef- 
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fective mean-field models to obtain the equation of state 
of hyperon rich matter up to densities relevant to the 
neutron star. The extra terms in the Lagrangian plays 
crucial role in determining the nature of an EoS, and 
is found to be in agreement with the Dirac-Bruckner- 
Hartree-Fock (DBHF) results at high density. The details 
of the hadronic model and the derivation of the equation 
of motion of the meson fields along with the energy den- 
sity and pressure for nuclear matter can be found in Ref. 
hj and we shall only briefly outline the formalism here. 
The field equations for the baryons and mesons are ob- 
tained from the Lagrangian density which is given by, 
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The scalar (psb) and the vector densities (p_b)j for a 
particular baryon species are given by, 



In the above lagrangian the covariant derivative I?^ is 
defined as Dp,= iguiB^p + igpBhBT°- p^^, where 
and fipi, are the field terms for the vector [lo) and the 
isovector meson (p) defined as follows. 



= dfi^u - d^cjfj,. (2) 
Rpu = di^p^^-dr^p^^+gpeabcplpt- (3) 

The above lagrangian (eqn. 1) embodies scalar meson 
interaction terms up to quartic order with K3 and K4 as 
the higher order scalar field constants for the cubic and 
the quartic terms respectively. Similarly for the vector 
(tj— meson) interactions, in addition to the usual second 
and fourth order mass terms, the scalar-vector cross cou- 
plings are also taken, where 771 and r]2 are the correspond- 
ing coupling constants. Similarly, the iso-scalar and iso- 
vector cross interaction term is also included (last line), 
where rjp is the constant that enters as a parameter in the 
model. Finally Co is the constant for the quartic vector 
interaction. Therefore, in addition to the scalar meson 
mass, and the three coupling constants namely, gcr,guj 
and gp, there are six other parameters that is tuned to 
reproduce the nuclear matter saturation properties. The 
influence of the inclusion of extra terms in the lagrangian 
on nuclear saturation properties can be seen in Table 1, 
where we have specified the models along with the corre- 
sponding saturation properties given in the lower panel. 

The meson fields equations for a, lo and p-meson, in 
the mean-field ansatz are given by 
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where Eg — y^k^ + m*^ is the effective energy and I3B 
is the isospin projection of B, the quantity fcs is the 
Fermi momentum for the baryon. The effective mass for 
the baryon species is then given by m*^ = mB — ga-BO'. 

The meson field equations for the cr, uj and p— mesons 
are then solved self-consistently at a fixed baryon den- 
sity to obtain the respective field strengths. The EoS 
for the /9— equilibrated for the hyperon rich matter is 
obtained with the requirements of conservation of total 
baryon number and charge neutrality condition given by, 



Y^QbPb+Y^QiPi^O, (9) 

B I 

where pB and pi are the baryon and the lepton (e,/i) num- 
ber densities with Qb and Qi as their respective electric 
charges. The lepton Fermi momenta are the positive real 
solutions of (/Cg -I- mD^^'^ = and (fc^ -I- m^)^/^ = Pp. = 

Me- 

The corresponding pressure and energy density of 
the charge neutral beta-equilibrated neutron star mat- 
ter(which includes the lowest lying octet of baryons) is 
then given by 
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The subscript '_B' in the above equations corresponds 
to the lowest lying octet of baryons and W and T spec- 
ifies the nucleon and leptons respectively. Here Pi and 
£1 are the contributions to the total pressure and energy 
density respectively from the leptonic counterparts. The 
corresponding model parameters that enters in our cal- 
culation of the equation of state at specified in Table I. 



TABLE L The model parameters from the relativistic mean- 
field theory. The lower panel display the nuclear matter sat- 
uration properties. 





G2 


Gl 


TMl* 


TMl 


NL3 


ma/mM 


0.554 


0.540 


0.545 


0.545 


0.541 


Qsn/AtX 


0.835 


0.785 


0.893 


0.798 


0.813 


Qvn/A-K 


1.016 


0.965 


1.192 


1.003 


1.024 


Qpn/Atv 


0.755 


0.698 


0.796 


0.737 


0.712 




3.247 


2.207 


2.513 


1.021 


1.465 


K4 


0.632 


-10.090 


8.970 


0.124 


-5.668 


Co 


2.642 


3.525 


3.600 


2.689 


0.0 




0.650 


0.071 


1.1 


0.0 


0.0 


m 


0.110 


-0.962 


0.1 


0.0 


0.0 


Tip 


0.390 


-0.272 


0.45 


0.0 


0.0 


a^MeV) 


-16.07 


-16.14 


-16.30 


-16.30 


-16.24 


Po(/m-^) 


0.153 


0.153 


0.145 


0.145 


0.148 


K (MeV) 


215 


215 


281.1 


281.1 


271.5 




0.664 


0.634 


0.634 


0.634 


0.595 


J (MeV) 


36.4 


38.5 


36.90 


36.90 


37.40 



From the table, one can clearly outline the differences 
in the nuclear matter saturation properties, which is at- 
tributed to the underlying interactions. For example, 
model Gl and G2, represents a class of model which en- 
tertains all the interactions described in the lagrangian 



and hence has the lowest value of incompressibility in 
comparison to other mean-field models given in Table I, 
i.e., the resulting EoS is softer. However, it is interesting 
to find that, although TMl* has the same attributes, yet 
the resulting incompressibility and all other saturation 
properties remains the same as TMl, where the former is 
a modification of the later with new couplings which stem 
from the modern effective field theory approach to rela- 
tivistic nuclear phenomenology. Among all the present 
models, NL3 has the lowest value of effective nucleon 
mass. It is to be noted that the K4 have a positive sign 
in TMl, TMl* and G2 unhke NL3 and Gl. The positive 
sign of K4 is required so that energy spectrum has a lower 
bound, else instabilities in calculations of the equation of 
state and of finite systems may occur [l^. In case of 
TMl* the inclusion of new couplings not only give the 
results close to DBHF at high density, it is also applica- 
ble to finite nuclei calculation for Z < % [lj[. It would 
be interesting to analyze the effect of these fundamental 
differences at high density and neutron star properties. 



III. STELLAR EQUATIONS 

The EoS of dense matter is the primary ingredient 
required to evaluate the global properties of the neu- 
tron stars. In order to account for the crustal part of 
the compact star, we include the BPS equation of state 
initially followed by the charge neutral, beta equi- 
librated hyperon rich matter, as described in the previ- 
ous section. However, in addition one needs to specify 
the hyperon coupling strengths . Presently, we assume 
that all the hyperons have the same couplings, which 
are in accordance with the quark sum rule approach. 
Accordingly, the scalar, vector and the iso-vector cou- 
pling strength for the hyperons are are expressed in terms 
of the nucleon coupling as — xh<j/xn<j = ■\/2/3, 

= XHui/xNw = a/s/S and Xp = XHp/xNp = \fW^- 
The equations for the structure of a relativistic spherical 
and static star composed of a perfect fluid were derived 
from Einstein's equations by Tolman, Oppenheimer and 
Volkoff [16], which are 



dP_ 
dr 



G{e + P][M + Anr^P] 
'7 (r - 2GM) 



dM 
dr 



47rr^e, 



(12) 



(13) 



with G as the gravitational constant and M(r) as the 
enclosed gravitational mass. We have used c — 1. Given 
an EoS, these equations can be integrated from the origin 
as an initial value problem for a given choice of central 
energy density, (sc)- The value of r (= R), where the 
pressure vanishes defines the surface of the star. We solve 
the above equations to study the structural properties of 
a static neutron star using the EoS for the electrically 
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charge neutral hyperonic dense matter [13, ■ In case 
of an evolved neutron star, the approximation of a zero 
temperature perfect fluid neutron star matter is reason- 
able (T^. Assuming a uniform rotation with static, axial 
symmetric space-time, the time translational invariant 
and axial-rotational invariant metric in spherical polar 
coordinates {t,r, 9^(j)) can be written as 

ds^ = -e^''dt^+e^°'{dr^+r^d0'^)+e^>^r^ sm^ 9 {dcjy^ujdtf , 

(14) 

where the metric functions v,a,P,U! depends only on r 
and 6. For a perfect fluid, the energy momentum tensor 
can be given by 

T^"" ^ Pg^"" + {P + e)u^'v!', (15) 
with the four-velocity 



= -==(1,0,0,17). (16) 

V 1 — w 

Here 

v = {VL-uj)rsuvee^-'', (17) 

is the proper velocity relative to an observer with zero 
angular velocity and is the angular velocity of the star 
measured from infinity. Now we can compute the Ein- 
stein field equations given by 



Rt^u - ■^gtj.uR = SttT^^ (18) 

(where R^i, is the Ricci tensor and and R is the scalar 
curvature and with c = G = 1). A maximum limit for 
the stable rotation of a star fi/^, is set by the onset of 
mass shedding from the equator of the star. General 
relativistic expression for this Keplerian frequency TJ/f ' 
can be obtained by using the extremal principle to the 
circular orbit of a point particle rotating at the equator 
of the star . We use the code written by Stergioulas 
[2H based on the Komatsu-Eriguchi-Hachisu method to 
construct uniformly rotating star models. 

IV. RESULTS AND DISCUSSION 

The extreme of matter density and high pressure pre- 
vailing in the core of neutron star pave way for exotics 
such as hyperons or quarks. In the present work, how- 
ever we are concerned only with the aspects and possi- 
ble implications of the presence of hyperons (lowest ly- 
ing octet of baryons) in the core of these stars. The 
ambiguity or lack of precise experimental data for the 
hyperon-hyperon or hyperon-meson interaction compels 



us to take certain assumption. Therefore, presently we 
take the hyperon-meson couplings in accordance with the 
quark sum rule approach as described in the previous sec- 
tion and is taken to be same for all the hyperon species. 
Although few hypernuclei experiments exist, but only for 
the case of A'^, which shows bound state of the particle 
at normal nuclear density, which comes out to be « - 
30 MeV. Theoretically, Hyperons are known to be stable 
and dominant species in the matter composition at high 
pressure and dense environment. Their appearance and 
abundance in the matter is known to have substantial 
contribution to the energy of the system rather than the 
pressure, which gives rise to a net softening effect on the 
EoS. Quantitatively, the amount of softening depends on 
factors, such as the underlying model attributes at satu- 
ration density and on the choice of coupling taken. How- 
ever, previous studies with hyperons in dense matter and 
neutron stars seems to agree on two important aspects, 
each of which has important implications to the neutron 
star structure and composition. One of which is the ap- 
pearance and abundance of hyperon species and the sec- 
ond one is the resulting deleptonization of dense matter. 
Leptons primarily gets used up to maintain the charge 
neutrality of the matter leading to the deleptonization, 
which does not happen in the case of nucleon matter com- 
position. To account for the crustal effect, we included 
the BPS EoS at subnuclear densities to the EoS of the 
hyperon rich neutron star matter. 

The resulting /?— equilibrated EoS for hyperon rich 
matter is shown in Fig. 1(A). From the plot, it can 
be seen that NL3 which has the lowest nucleon effec- 
tive mass predicts stiff EoS than the rest of the models. 
Although not much of a difference is seen in the energy 
per nucleon as well as saturation density for the models 
, there is substantial difference in the incompressibility 
and nucleon effective mass on account of its sensitivity 
to the EoS. However, beyond the normal nuclear mat- 
ter density, the EoS is sensitive to the appearance and 
abundance of the hyperon species. Fig. 1(B) shows the 
density at which the first two members of the lowest lying 
octet of barons (A°(1116) and S~(1193)) start appearing 
in dense matter. The plot reveals that although is 
w 77 MeV massive than A", still they are the first ones 
to appear in dense matter owing to their negative isospin. 

seems to start appearing at nearly (1.5— 1.7)po closely 
followed by A*' at (1.7 — 2.2) po- The resulting softening 
effect on the EoS is reflected on the EoS for the models 
at £ « (200 - 250)MeT^/m~^. The successive appear- 
ance of other hyperon species adds to the softening of the 
equation of state thereafter. 

With the five EoS from the relativistic mean-field mod- 
els as input, we now calculate the corresponding stellar 
properties for static case. Fig. 2(A), shows the mass 
of the neutron star obtained as a function of the central 
density of the star in the static limit. The filled circles 
denote the point of maximum mass for the corresponding 
models, presently under investigation. The quantitative 
impact of the onset of hyperons on the EoS and hence on 





FIG. 1: (Color online) (A): Equation of state (e vs. p) of hy- 
peron rich neutron star matter for the five models from the 
relativistic mean- field theory. (B): Normalised baryon densi- 
ties (ps/po) at which E~ and A" appears in dense matter, for 
the corresponding models. 



FIG. 2: (Color online) For static case (A): Mass of the neu- 
tron star sequence plotted as a function of the central density 
of the star. (B): Mass of the neutron star as a function of 
radius is shown along with the constraints from neutron star 
observation and causality condition. 



the neutron star mass can be known if we compare the 
densities at which they appear. For example, 'E~' and 
'A*^' starts appearing nearly at similar density in case of 
TMl* and G2 and hence the resulting EoS follows simi- 
lar patters although G2 predicts the softest EoS among 
the five models considered here. Because of the similar 
behavior and softening at high density, they also result 
in star with high central density (Ec)- For example, in 
case of TMl*, Sc « 9po and for G2, it is « lOpo- For 
the present models, the central density of the star at 
maximum mass lies in the range (4 — 10) po- The cor- 
responding sequence of neutron stars in the M-R plane 
for static case is shown in Fig. 2(B) along with avail- 
able constraints. The maximum mass of the static star 
lies in the range (1.29— 1.79) Mq and the corresponding 
compactness ratio (M/R) is in the range (0.11 — 0.14). It 
is interesting to see that, the results compare very well 
with the observation of massive stars (M « 2Mq) [2^.[23j 
from observation of QPO's (Quasi Periodic Oscillations) 
from X-ray emmisions. However, there are also observa- 
tion of least massive stars with mass 1.18 ± O.O2M0 [23 | 



from the binary pulsar J1756— 2251. These lower and the 
upper bound on mass from observations have raised in- 
teresting implications on their structure and composition 
[2^ . However, the canonical value of M = 1.44M0 still 
remains the largest precisely measured mass of the Neu- 
tron star P5'i?1913-|-16 [2^. Under these circumstances, 
precise radius measurement holds the key to understand 
the structural properties of these stars. From the obser- 
vational point of view, there are large uncertainties in 
determination of the radius of the star [23, [H, [2^ which 
is primarily because of our lack of knowledge of the com- 
position of the star atmosphere, large distance and also 
due to the presence of high magnetic fields. However, it 
is worth to recall that although there are large obser- 
vational errors in the mass-radius determination, yet in 
case of VelaX — 1 the lower mass limit « (1.6 — 1.7)AfQ 
is at least mildly constrained by geometry [30] ■ Theoret- 
ically the mass of these compact stars are controlled by 
the nature of the EoS (stiffness/ softness) at an excess 
of 5poi whereas the determination of radius is controlled 
by the EoS in the vicinity of po, the nuclear matter sat- 



uration density. The region excluded by the causahty 
constraint i.e., ^Jdp/de < 1 leads to i? > 2.9GM/c^ is 
shown in Fig. 2(B) along with the redshift measurement 
from active X-ray burster EXO 0748-676, which implies 
z ~ 0.35 [1^. Although this value may not be universally 
accepted, yet EXO 0748-676 seems to be a promising can- 
didate as it can provide stringent constraint in the M-R 
plane. The experimental constraints are extracted from 
observations of kilohertz quasi-periodic brightness oscil- 
lations in the Low Mass X-Binaries 4U 0614-1-09 and 4U 
1636-536. Here we find that G2 results in low mass with 
highest central density owing to the soft EoS due to the 
underlying cross couplings included in the interactions. 
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FIG. 3: (Color online) For star rotating at kepler velocity (A): 
Mass of the neutron star sequence plotted as a function of the 
central density of the star. (B): Mass as a function of radius 
for the corresponding sequence. The filled circles denotes the 
values at maximum mass. 

We now calculate the global properties of the neutron 
star undergoing keplerian rotation. Compared to static 
case, the effect of rotation is to increase the equatorial 
radius of the star and also to increase the mass that can 
be sustained at a given central energy density. For star 
rotating at kepler velocity we show the obtained mass of 
the neutron star sequence plotted as a function of the 



FIG. 4: (Color online) For star rotating at kepler velocity 
(A): Variation of the Moment of Inertia (10*^(?cm^)' of 
the neutron star sequence as a function of the radius of the 
star. (B): The corresponding gravitational energy of the star 
as a function of mass is shown. (C) The angular momentum 
of the star is shown as a function of mass for the five model 
under consideration. The filled circles denotes the values at 
maximum mass. 



central density of the star in Fig. 3(A). We find that 
the mass of the neutron star increases for all the models, 
however the resulting central density of the star drops. 
Quantitively, we find an increase of « (15 — 20) % in the 
stellar mass and the decrease in the central density is « 
(5 — 17) % which corresponds to « (0.2 — 1.6)po. At the 
mass-shedding limit, the maximum mass of the star is 
found in the range (1.49 - 2.18) Mq. Fig. 3(B) shows 
the corresponding models of neutron stars at = ilfe 
(mass-shedding limit) in the M-R plane. The filled cir- 
cles shows the point of maximum mass for the models. 
In comparison to the static star, we find an increase of 
~ (38 — 43)% in the equatorial radius at maximum mass. 
The increase in the mass or equatorial radius are EoS 
dependent and in turn sensitive to the underlying model 
attributes at po- For example, in the framework of an 
effective chiral model, the neutron star properties in the 
kepler limit was studied and there an increase of nearly 
50% was found in the star radius [3l|- This was a con- 
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FIG. 5: (Color online) (A): The kepler velocity of the star 
is plotted as a function of mass of the star. (B): The kepler 
velocity of the star as a function of the resulting deformation 
parameter {Rp/Re). The filled circles denotes the values at 
maximum mass. The rotational frequency of the two fastest 
known pulsars, rotating at 1122Hz and ll&Hz are also indi- 
cated (dashed horizontal lines). 

sequence of the underlying softness of the EoS. Similar 
correlations can be found in other related properties of 
the rapid rotors, which we discuss next. 

Fig. 4(A) display the moment of inertia of the sequence 
of rapid rotors as a function of the radius of the star. Mo- 
ment of inertia is a sensitive quantity for rapid rotors in 
addition to the fact that its accurate determination from 
observations would imply a higher degree of accuracy in 
the determination of neutron star radius, which still re- 
mains a major hurdle in our present understanding of 
these stars. Dimensionally, / c>c MR^ and they play 
crucial role in the models of radio pulsar. Independent 
of the rotation, i.e., slow or fast, the relation of moment 
of inertia to the matter distribution within the star is 
complicated. For the models presently under investiga- 
tion, the moment of inertia of the compact stars lies in 
the range 1^^ = (1.5 — 3.9), where G2 results in lower 
value and NL3 has the highest value and has the largest 
compactness parameter (M/R) among them. The general 
relativistic limit for the compactness parameter assuming 



a uniform density star with the causal equation of state 
i.e., P = e gives M/R < 4/9 [S^l- The compactness ratio 
for the present models lies far below the aforesaid limit 
and is found in the range M/R = (0.10 — 0.12). From 
the observational point of view, recently discovered rela- 
tivistic pulsar PSRJ0737 - 3039 [H could be the first 
one with moment of inertia measured precisely. However, 
the estimates from crab pulsar [33| puts 745 > 1.61 for 
Mneb = 2.OM0, which is considered to be the conserva- 
tive estimate, and I45 > 3.04 for Mmb = 4.6M0 (latest 
estimate). Fig. 4(B) and Fig. 4(C) shows the variation 
of the gravitational/ rotational energy of the star rotat- 
ing at kepler velocity as a function of stellar mass and 
the corresponding angular momentum of the rapid ro- 
tors respectively. As anticipated NL3 stands apart from 
the rest of the models, which reflects the stiffness of the 
corresponding EoS. 

Fig. 5(A) shows the variation of kepler velocity of the 
star as a function of star mass. Initially, we find that 
the rotation of the star increases linearly with mass for 
all the cases, but after the star attains a maximum mass 
of « (1.4 — 1.5)M0, they undergo much rapid rotation 
but the increase in mass is not so substantial with an ex- 
ception of NL3. In the inset on y axis, the corresponding 
rotational period is shown, which indicates that the stars 
composed of hyperon rich matter can rotate much faster 
to be observed as a millisecond pulsar. For all the mod- 
els presently studied, the neutron star seems to attain 
very high velocity in the range fife = (6500 — 7500)s~^ 
at maximum mass, which corresponds to a kepler period 
of Pfc = (0.966 - 0.837). It is interesting to flnd that 
the stars corresponding to the softest two prescription 
i.e., Gl & G2 attains minimal and maximal rotation re- 
spectively. Recently there has been considerable debate 
over the possibility of finding a millisecond pulsar and 
its implication on the star structure. In the Fig. 5(A), 
we plot the two fastest known pulsars till date and com- 
pare it with the results obtained in our calculations. The 
observation of rapid rotors such as XTE J1739-285 with 
spin frequency of 1122 Hz is the most rapidly rotating 
neutron star yet discovered [sl]. Prior to the discovery 
of 1122 Hz, the discovery of a radio pulsar spinning at 
716 Hz suggests that the true maximum spin rate can be 
even higher j36i] . However, there are few more observa- 
tions of pulsar period rotating in the vicinity, namely the 
PSRB1937 + 2Ui^ = 633 Hz) ^ and PSRB1957 + 20 
{ly — 621 Hz) [38|. It is noteworthy that the minimum 
observed pulsar period is still 1.56 ms and it is often ar- 
gued that various mechanisms such as the r-mode insta- 
bility and gravitational radiation may help in the slow 
down of the pulsar velocity and hence the millisecond 
pulsars are an exception. Such high rotations would in- 
variably throw interesting possibilities pertaining to the 
structure and composition of a compact star, which we 
discuss next. 

Rapid rotation can cause substantial deformation in 
the stellar structure as the matter is thrown in the equa- 
torial plane, thereby resulting in flattening effect at the 
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FIG. 6: (Color online) The stellar structure for the static 
(n = 0; left one) and the corresponding ones rotating at kepler 
velocity (f2 = O^.; right one) is shown for the models studied 
in the coordinates of polar radius (y axis) and the equatorial 
radius (x axis) in km. 



TABLE II: Comparison of the global properties of the neutron 
star in the static limit i.e., = (upper row) and the stars 
rotating with kepler velocity (n = ilk) at maximum mass 
(lower row) for the corresponding EoS as tabulated. 



EoS 


M 


Re 


Sc/so 


Mb 


I45 


/ 


p 1 ^ 




{Mq ) 


(Km) 




(Mo) 


iW°gcm^) (Hz) 




NL3 


1.79 


13.20 


4.64 


2.00 










2.18 


18.23 


4.40 


2.44 


3.92 


1095 


0.56 


TMl 


1.51 


12.36 


6.17 


1.65 










1.80 


17.06 


6.00 


1.97 


2.49 


1100 


0.58 


TMl* 


1.42 


11.30 


8.99 


1.55 










1.68 


16.50 


6.59 


1.83 


2.11 


1117 


0.58 


Gl 


1.33 


12.45 


5.06 


1.44 










1.58 


17.33 


5.00 


1.70 


2.05 


1009 


0.59 


G2 


1.29 


10.64 


10.15 


1.40 










1.48 


15.25 


8.42 


1.61 


1.47 


1183 


0.59 



poles. On similar grounds, the deerease in the central 
density of the star can be understood in case of fast 
rotors. From the plot (Fig. 5(B)), we find that NL3, 
which results in the stiffest EoS among the five models 
studied here, undergoes maximum deformation and G2 
is the least deformed star. However, for all the mod- 
els the star seems to undergo maximal deformation at 
Vtk w 4700s"\ (/ » 716Hz) thereafter the star ro- 
tates at much higher frequency without any further de- 
formation (an increase in velocity of w 50%). The shape 
of the star for in the static limit and the corresponding 
star rotating with kepler velocity at maximum mass ob- 
tained, is shown in Fig. 6, for the five models in the 
polar-equatorial plane. We find that G2 sustains max- 
imal keplerian velocity and the Gl is credited with the 
least, although the resulting deformation for the two re- 
mains the same. However, at maximum mass the star 
rotates with nearly equal frequency in case of TMl and 
TMl*. As expected, the stiffness of the EoS showed up 
in case of NL3, which results in maximal deformation 
and G2 deforms the least. For a better comparison, the 
global properties of the neutron star with hyperon core, 
corresponding to the static and the kepler limit for the 
models under investigation are tabulated in Table II. 



V. CONCLUSIONS 

We study the effect of nuclear interaction terms on the 
equation of state of hyperon rich matter and the resulting 
consequence on the neutron star properties in both the 
static and the kepler limit. The effect of rapid rotation on 
the stellar structure is investigated with different models 
and the effect of rotation on the neutron star structure is 
lineated. Rotating neutron stars, can potentially not only 
constrain the EoS of dense matter, but is also interest- 
ing because they can be the source of gravitational wave 
radiation, the mechanism which is believed to take away 
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the star angular momentum thereby resulting in the slow 
down of the stars. For the purpose, we took five different 
models from effective mean-field theory with different in- 
herent features, such as the model interactions and the 
nuclear saturation properties to analyze their effect on 
the equation of state of dense matter. The inclusion of 
extra terms in the lagrangian contributes to the softening 
of the EoS and the matter becomes more compressible. 
From the result of the properties of the star, we found 
that a soft EoS (one with lower value of incompressibility) 
results in star with small maximum mass, small radius, 
and large rotation rate and vice versa. In comparison 
to the static stars, the rotating stars gain more mass. 
We found an increase of « (15-20)% in the maximum 
mass, with a corresponding decrease of « (5-17) % in 
the central density of the star rotating with kepler ve- 
locity. The models can also explain fast rotors, typically 
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